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PREFACE  | 

This  investigation  was  sponsored  by  the  Defense  Nuclear  Agency 
under  Subtask  Y99QAXSB211,  "Penetration,"  Work  Unit  20,  "Penetration 
Support."  This  study  was  conducted  by  personnel  of  the  Structures 
Laboratory  (SL),  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES), 
during  October  1978  through  May  1979,  under  the  general  supervision  of 

i 

Mr.  Bryant  Mather,  acting  Chief,  SL,  and  Dr.  J.  0.  Jackson,  Jr.,  Chief, 

Geomechanics  Division  (GD),  SL.  Mr.  R.  S.  Bernard  formulated  the  < 

! 

theory,  and  Mr.  D.  C.  Creighton  developed  and  implemented  the  computer  j 

analysis,  both  under  the  technical  guidance  and  direction  of  Dr.  B.  Rohani, 

GD.  Messrs.  Bernard  and  Creighton  prepared  this  report.  4 

COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE,  were  Command¬ 
ers  and  Directors  of  WES  during  the  period  of  research  and  report  prepara- 
tion.  Mr.  F.  R.  Brown  was  Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  cf  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (Si)  units  as  follows: 


Multiply  By 

feet  0.301*8 

inches  0.025** 

feet  per  second  0.301*8 

pounds  (force)  per  foot  lU.5939 

pounds  (force)  per  square  inch  689U.757 
pounds  (mass)  0.1*535921* 

pounds  (mass)  per  cubic  foot  16.0181*6 

pounds  (mass)  per  inch  17*85797 

pounds  (mass )-square  inches  0.0002926 


_ To  Obtain _ 

metres 

metres 

metres  per  second 
newtons  per  metre 
pascals 
kilograms 

kilograms  per  cubic  metre 
kilograms  per  metre 
kilograms-square  metres 
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PROJECTILE  PENETRATION  IN  SOIL  AND  ROCK: 
ANALYSIS  FOR  NON-NORMAL  IMPACT 


CHAPTER  1 
INTRODUCTION 

1.1  BACKGROUND 

In  recent  years,  due  to  growth  of  the  experimental  data  base, 
development  of  empirical  formulas,  and  implementation  of  new  theoretical 
and  computational  techniques,  it  has  become  possible  to  predict  the 
approximate  motion  of  stable  earth  penetrators  after  normal  impact  in 
many  types  of  geological  targets  (Reference  1-1).  Results  of  some  non¬ 
normal  (yawed/oblique)  impact  tests  have  also  been  predicted  success¬ 
fully,  although  the  data  base  is  small  compared  to  that  for  normal  impact 
(References  5  through  7). 

Previous  analyses  by  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  of  non-normal  impact  have  used  the  Cavity  Expansion  Theory 
(CET)  for  rock,  and  an  empirical  formulation  similar  to  Young’s  equation 
for  soil  (References  5  and  8).  These  analyses  were  considered  prelimi¬ 
nary  when  they  were  first  developed,  and  it  was  hoped  that  more  credible 
analyses  would  arise  in  the  future. 

1.2  PURPOSE  AND  SCOPE 

The  present  investigation  concerns  the  development  of  a  generalised 
analysis  for  oblique/yawed  impact  and  penetration  in  soil  and  rock  (or 
rock-like)  targets,  based  on  the  modification  of  previous  analyses  and 
the  interpretation  of  recent  test  data.  The  objectives  are  (a)  to 
obtain  external  load-time  histories  sufficiently  accurate  for  structural 
response  calculations  and  (b)  to  calculate  postimpact  trajectories  well 
enough  to  predict  the  terradynami c  performance  ami  stability  of  a  given 
projectile . 

The  two-dimensional  (2-D)  penetration  theory,  including  equation  of 
motion,  stress  distribution,  and  free-surface  and  wake  separation- 
reattachment  effects,  was  developed.  The  penetration  theory  was  computer 

T 
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coded  into  PENC02D  to  make  comparison  calculations  with  experimental 
data  from  both  reverse  ballistic  and  conventional  penetration  tests. 
Using  PENC02D,  the  effects  of  initial  impact  conditions,  projectile 
geometry,  and  soil  penetrability  were  examined  in  a  brief  parameter 
study. 
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PENETRATION  ANALYSIS  FOR  MOTION  IN  TWO  DIMENSIONS 


2.1  PROJECTILE  MOTION  IN  TWO  DIMENSIONS 

Consider  a  rigid,  deeply  buried  projectile  (no  free-surface  effects), 
whose  motion  lies  in  the  XZ-plane,*  as  shown  in  Figure  2.1.  The  rota¬ 
tion  is  described  by  the  angular  velocity  0  ,  and  the  translation  by 

•  • 

X  .  j.  Z  ,  the  X-  and  Z-velocity  components  of  the  center  of  gravity 

(CG).  For  analysis,  it  is  convenient  to  express  the  CG  velocity  in 

terms  of  its  lateral  and  axial  (x-  and  z-)  components  V  and  V  , 
b  x  z 

respectively: 


V  =  X  cos  0  -  Z  sin  0 
x 


(2.1) 


V  *  X  sin  0  +  Z  cos  0 
z 


For  purely  axial  motion  (V^  »  0  ,  0=0),  the  stress  distribution  on 

the  projectile  is  symmetric,  and  the  net  force  is  in  the  axial  direc¬ 
tion,  producing  no  pitch  (turning)  moment.  The  introduction  of  rotation 
(0  ^  0)  or  lateral  motion  (V^  #  0)  destroys  this  symmetry,  and  it  is 
necessary  to  specify  the  asymmetric  stress  distribution  in  order  to 
calculate  the  resulting  force  and  moment. 

Suppose  that  the  compressive  resisting  stress  a  is  always  normal 
to  the  projectile  surface  (no  tangential  stresses).  On  any  surface  area 
element  dA  ,  the  x-  and  z-components  of  o  are,  respectively. 


2  The  XZ-coordinate  system  is  fixed  in  the  target  (Figure  2.1). 

A  dot  above  any  quantity  indicates  differentiation  with  respect  to 
^  time. 

For  convenience,  symbols  and  abbreviations  used  in  this  report  are 
listed  and  defined  in  the  Notation,  Appendix  C. 

The  xz-coordinate  system  is  fixed  on  the  projectile  axis,  with  its 
origin  at  the  CG  (Figure  2.1).  Equations  2.1  and  2.2  represent 
the  projections  of  the  absolute  velocity  vector  along  the  x-  and 
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( 


a  •  a  cos  n  cos  * 
x 


(2.3) 


a  ■  a  ain  n 
2 


(2.M 


where 

p  *  asimuthal  angle  (Figure  2.2) 

tan  n  ■  olope  of  projectile  surface  at  a  given  point,  with  respect 
to  axis  of  symmetry  (Figures  2.1  ami  2.3) 

The  lateral  and  axial  force  components  (Figure  2.2)  acting  on  surface 
element  dA  are,  res  'cctlvely , 


dF  ■  ~o  dA 
x  x 


(2.5) 


dF  «  -o  dA 
s  a 


where  dA  in  the  differential  element,  of  surface  area.  The  lateral 
and  axial  net  force  components  acting  on  the  entire  projectile  are  then 


F  -  o  dA 
x  x 


(2.7) 


F  *  -/  o  dA 
a  s 


where  the  Integral  is  evaluated  over  the  entire  projectile  surface. 
The  translational  (Oil)  equations  of  motion  are 


M  X  *  F  cos  0  >  F  sin  P 
x  a 


M  *  -F  sin  0  >  F  cos  0 

x  7. 


(2. 10) 


where  M  Is  t  lie  projectile  mass  and  X  mid  are  the  components  of 
acceleration  in  the  X-  and  '.'.-direct  ions,  respectively.  The  projectile 
is  assumed  to  be  rigid,  so  the  rotational  equation  of  mot  Ion  Is 


IP  -  /  7.  dF  -  /  x  dF 

X  5 


(2.11) 
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where  I  in  the  transverse  mean  moment,  of  inertia  about  the  Cl! ,  0  la 

the  angular  acceleration,  ami  x  ami  ?,  are  the  lateral  and  axial 
distances,  respect Ively ,  from  the  CO  (Figure  2.1).  The  right-hand  side 
of  Equation  2.11  Is  the  total  moment  exerted  on  the  projectile  about 
its  CG,  and  the  Integrals  are  evaluated  over  the  entire  area  of  the 
pro,)  ec  1 1 1  e  s  u  r  f ac  e . 

2 . 2  CTHECo  nUlTHIlUlTION  FOR 

ROCK  AND  CONCRETE 

For  penetration  Into  concrete  and  Intact  rock  (Appendix  A),  the 
net  force  on  the  projectile  can  be  calculated  approximately  from  the 
following  normal  stress  distribution: 


where 

Y  »  uncon fined  compressive  strength  of  Intact  target  material  . 

C  »  target  mass  density 

v  «  absolute  local  velocity  <'f  a  given  point  on  projectile  surface 
(Figure  2.  j) 

v  ’  outward  normal  component  of  v  l Figure  2.  0 

IVflnml  In  t  arwn  of  ,  0  %  ami  a  ,  thr  qunnt.it.  iru  v  anil 

v  nro 
n 


-  (v 


(2.1  0 


If  necessary,  minor  correct  ions  !  o  Y  can  be  made  for  the  quality 
of  the  rock  or  for  t  lie  maximum  aggregate  sire  In  the  oonorete 
(Appendix  A).  Any  set  of  units  can  be  used  in  Equation  2.12  but 
they  must  be  dimensionally  compatible  (e. g. ,  s lugs ,  feet, 
seconds! . 


1 1 
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<2. 14) 


v  *  V  sin  n  ♦  (V  +z0)  cos  n  cob  6 
n  z  x 


It  ia  aa a  tuned  that  a  given  surface  element  must  have  a  net  motion  Into 
the  target  material  to  produce  a  reoisting  stress;  hence,  o  vanishes 
when  vn  0  in  Equation  2.12. 

2.3  STRESS  DISTRIBUTION  FOR  SOIL 

For  penetration  into  soil  (Appendix  B),  the  net  force  on  the 
projectile  can  be  calculated  approximately  from  the  following  normal 
stress  distribution: ^ 


(2.15) 


I 

I 

i 


where 

U  ■  1.96  *  10^  lb/ft,  constant^ 

/  o  0  7 

6  *  6.25  x  10  lb  -sec  /ft  ,  constant 
Y  ■  1.56  *  10^  lb/ft^,  constant 

S  *  Young's  S-number  (Appendix  B),  soil  penetrability  index 

r  ■  local  cylindrical  rndius  at  a  given  point  on  the  projectile 
”  surface 

Z  -  vertical  depth  of  a  given  point  on  the  projectile  surface, 
measured  from  the  target  surface 

The  expressions  for  v  and  vn  are  still  given  by  Equations  2.13  and 
2.14,  respectively;  and  the  requirement  that  v  >0  for  0  >  0  applies 
in  soil  as  well  as  in  rock  or  concrete. 

The  quantities  u  ,  0  ,  and  y  are  constants,  independent  of  soil 


Any  set  of  dimensionally  compatible  unit3  can  be  used  in  Equation  2.15, 

■j  as  long  as  all  quantities  are  converted  to  those  units.  j 

A  table  of  factors  for  converting  U.  S.  customary  units  of  measurement  > 

to  metric  (SI)  units  is  presented  on  page  h.  i 
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type  and  projectile  characteristics.  Thus,  the  only  parameter  charac¬ 
terizing  the  soil  is  Young's  empirical  S-number,  which  can  be  obtained 
for  a  given  target  from  previous  penetration  data  or  from  correlation 
with  geological  descriptions  (Table  B.l). 

2.U  WAICE  SEPARATION  AND  REATTACHMENT 

Flash  X-ray  photographs  and  two-dimensional  (2-D)  finite-difference 
calculations  have  shown  that  loss  of  contact  between  target  and  projec¬ 
tile  occurs  somewhere  on  the  nose,  while  reattachment  may  or  may  not 
occur  on  the  aftbody.  The  current  analysis  is  not  sophisticated  enough 
to  predict  the  actual  onset  of  separation.  Nevertheless,  it  is  neces¬ 
sary  to  include  a  simple  model  for  the  kinematics  of  wake  separation, 
especially  when  2-D  projectile  motion  in  soil  is  investigated. 

It  is  assumed  that  there  is  a  minimum  angle  of  approach,  or  wake 
separation  angle,  required  between  the  target  and  the  projectile 

contact  surface  (Figure  2. h)  in  order  for  contact  to  be  maintained. 

For  a  given  surface  element,  the  local  angle  of  approach  4>  is  deter- 

g 

mined  by  the  instantaneous  orientation  of  the  projectile  such  that 

V 

sin  $  *  —■  (2.16) 

where  V  is  the  CG  velocity  and  Vn  is  its  outward  component  normal 
to  the  projectile  surface,  given  by 

V  =  V  sin  n  +  V  cos  n  cos  ip  (2.17) 

n  z  x 

Thus,  separation  occurs  whenever  $  <  <p  .  ,  i.e.,  whenever  the  local 

-  min  q 

angle  of  approach  is  less  than  the  wake  separation  angle. 

Separation  is  only  part  of  the  problem;  reattachment  may  or  may 

8  ■ 

The  rotational  velocity  9  is  not  considered  in  formulating  the 
criterion  for  wake  separation.  This  quantity  does,  however,  play 
q  a  strong  role  in  wake  reattachment. 

For  the  case  of  no  yaw  (no  angle  of  attack),  Eqxiation  2.16  reduces 
to  sin  $  *  sin  n  . 
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not  occur,  depending  on  the  relative  motion  of  the  wake  and  the  aftbody 
downstream  of  the  separation  point.  It  is  expected  that  the  wake  separa¬ 
tion  angle  will  generally  he  small  (<1C  degrees),  in  which  case 

the  angle  of  attack  a  required  to  maintain  contact  will  also  he  small. 
In  any  case,  reattachraent  will  he  analyzed  based  on  three  assumptions: 

1.  Constant  axial  velocity  (V^) 

2.  Constant  angular  velocity  (6) 

3.  Constant  angle  of  attack  (a) 

Obviously,  these  three  quantities  change  as  the  projectile  penetrates 
into  the  target.  Nevertheless,  if  the  quantities  do  not  vary  rapidly, 
the  quasistatic  approximation  is  adequate. 

The  cavity  formed  by  the  wake  is  a  gently  curved  cylindrical  tube 
of  increasing  radius  (with  time),  and  the  projectile  rotates  inside 
this  tube  a3  it  travels  forward.  The  curved  axis  of  the  cavity  is 
fixed  in  the  target,  and  aftbody  reattachment  occurs  wherever  the 
projectile  rotates  into  the  cavity  wall  (Figure  2.5). 

Disregarding  any  stresses  that  may  oppose  the  radial  expansion  of 
the  cavity,  conservation  of  mass  and  incompressibility  in  the  target 
requires  that 

r  r  ■  constant  (2.18) 

c  c 

where  r  is  the  local  cylindrical  wake-cavity  radiu3.  Evaluating 
c 

Equation  2.18  at  the  separation  point,  it  then  follows  that 


r  r  *  r  V  tan 
c  c  o  z 


(2.19) 


where  rQ  is  the  projectile  radius  at  the  separation  point.  For  small 

values  of  $  .  ,  r  is  approximately  the  radius  corresponding  to 

min  o 

n  *  (Figure  2.1*).  Furthermore,  for  straight-  and  tapered-aftbody 

projectiles,  rQ  is  approximately  the  radius  at  the  base  of  the  nose. 

Denoting  axial  distance  aft  of  the  nose  tip  by  c  ,  time  deriva¬ 
tives  can  be  related  to  V  and  C  by  the  transformation 

z 


(2.20) 


<L-,v  i- 

dt  z  dC 


and  Equation  2.19  reduces  to 


drc 

r  -r—  =  r  tan  $  . 
c  d£  o  min 


(2.21) 


for  which  the  solution  is 


s  Vr o2  +  2ro  ?  -  ?o)  tan  Vin  ’  C  i  ?c 


(2.22) 


where  _ is  the  value  of  c  at  the  separation  point.  For  ?  <  5  , 

the  cavity  radius  is  the  same  as  the  projectile  radius,  since  no  loss 
of  contact  has  occurred. 

Using  Equation  2.20  again,  the  differential  equation  for  the 
lateral  displacement  5  of  a  point  on  the  projectile  axis,  relative 
to  the  cavity  centerline  (Figure  2.5),  is 


V  Si-  r  .  e  0 

dt  z  d C  ^  S 


(2.23) 


Treating  V  and  9  as  constants,  the  solution  for  5  is 
z 


(2.2U) 


The  initial  displacement  5q  =  5  sin  a  is  the  displacement  due  to  the 
angle  of  attack  a  (treated  here  as  a  constant ) ,^  so  Equation  2.24  is 
replaced  hy 
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It  is  assumed  that  small  attack  angles  do  not  significantly  change 
the  cavity  geometry.  Thus,  if  9=0,  the  cavity  centerline  will 
he  straight;  hut  if  in  addition  a  f  0  ,  the  projectile  axis  will 
he  skewed  relative  to  the  cavity,  producing  displacements  5  = 

C  sin  a  between  the  projectile  axis  and  the  cavity  centerline. 
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(2.25) 


where 


C  sin  a 


a  ■ 


(2.26) 


At  any  point  on  the  projectile  surface  aft  of  the  separation  point,  the 
criterion  for  no  contact  is1* 


2  (c  "  co)2  v~  *  5  sin  a  '  rp  COB  *  <  ^rc2  "  rp2  8in?  *  ’  2  -  *  -  21 

■  2  (;  -  <o)2  r-'slna<  / 


(2.27) 


rp  cos  * 


2  2.2.  ir  .  » 

r  -r  a  in  *  ,  -:r  <  ♦  <  :r 

c  p  2  —  —  2 


Hence,  if  ^  <_  and  Equation  2.27  are  satisfied,  there  is  no  con¬ 

tact,  and  o*0. 

The  foregoing  analysis  is  only  a  rough  method  of  accounting  for 
wake  separation  and  re attachment  using  a  single  input  parameter,  the 
wake  separation  angle  •  Nevertheless ,  this  analysis  does  make 

it  possible  to  assess  the  relative  effects  of  separation  and  reattach¬ 
ment  on  projectile  stability,  as  will  be  shown  in  Chapter  3. 


2.5  FREE- SURFACE  EFFECT 

In  brittle  materials,  such  as  rock  and  concrete,  a  crater  is  always 
formed  during  (or  Just  after)  impact.  Postirapact  inspections  in  rock 
(Reference  9)  indicate  that  these  craters  may  be  as  much  as  10  projec¬ 
tile  diameters  across  at  the  target  surface.  The  crater  width  decreases 


Figure  2.6  presents  a  rear  axial  viev  of  the  projectile  and  cavity 
cross  sections  at  an  arbitrary  location  aft  of  the  separation  point. 
The  figure  is  valid  only  for  small  displacements,  i.e.  6  «  t  . 


16 


1 

rapidly  with  depth,  however,  approaching  a  value  equal  to  (or  slightly 

less  than)  the  projectile  diameter  after  a  few  calibers  of  penetration. 

The  same  phenomenon  occurs  in  soil,  though  to  a  lesser  extent  than  in 

rock,  since  soils  are  more  ductile  than  brittle  in  behavior. 

For  normal  impact  or  for  near-normal  impact  with  slight  yaw,  it  is 

not  necessary  to  account  for  the  cratering  phenomenon  (i.e.,  free- 

surface  effect)  in  calculating  the  loads  on  the  projectile  with  the 

current  analysis.  On  the  other  hand,  for  oblique  impact  with  no  initial 

yaw,  the  free-surface  effect  is  probably  one  of  the  dominant  mechanisms 

creating  an  unbalanced  lateral  force  on  the  projectile,  especially  in 

rock  and  concrete.  Thus,  for  analyzing  non-normal  impact  in  general, 

the  presence  of  the  free  surface  should  be  acknowledged,  even  if  only 

by  a  crude  representation  of  its  effect. 

The  influence  of  the  free  surface  is  modeled  herein  with  an  on-off 

stress  criterion,  governed  by  the  location  of  a  given  projectile  surface 

element  dA  with  respect  to  the  free  surface  (Figure  2.7).  Whenever 

the  radial  distance  from  the  projectile  axis  (through  JA)  to  the  free 

12 

surface  is  less  than  some  prescribed  value  rg  ,  then  the  stress  on 
dA  is  set  equal  to  zero.  Quantitatively,  this  condition  is  expressed 
by  setting  a  0  whenever 

ZCG  ♦  z  cos  9  <  rg  sin  9  cos  <|>  (2.28) 

where  Z„„  is  the  vertical  distance  from  the  free  surface  to  the  CG. 

CG 

Assuming  rg  to  be  directly  proportional  to  the  local  projectile  radius 
Tp  ,  Equation  2.28  is  replaced  by  the  stipulation  that  a  =  0  whenever 

Z_„  +  z  cos  9  <  kr  sin  9  cos  <i>  (2.29) 

CU  p 

where  the  free-surface  parameter  is 

12 

r3  is  the  maximum  perpendicular  distance  from  projectile  axis, 
through  dA  ,  to  target  surface  for  which  stress  relief  due  to  the 
free  surface  can  occur. 
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(2.30) 


r 

k  ■  —  ■  constant 

rp 

This  model  for  the  free-surface  effect  produces  a  stress  distribu¬ 
tion  that  may  be  partially  or  completely  turned  off  on  the  top  side 
(-»/2  <  $  <  n/2)  of  the  projectile,  depending  on  the  geometry,  depth, 
and  orientation  vith  respect  to  the  target  surface.  The  deeper  the 
penetration  and  the  more  vertical  the  orientation,  the  smaller  the 
influence  of  the  free  surface  on  the  lateral  loads. 


Q  oc, 


\  /wafeble 


Figure  2.1  Projectile  motion  in  two  dimensions. 


Figure  2.2  Three-dimensional  view  of  projectile 


Figure  2.5  Separation  and  reattachment  on 
a  rotating  projectile. 


/ 


WAKC-CAVITV  WALL 


\ 

\ 


I 

/ 


NOTE.  THIS  DRAWING  VALID  ONLT  FOR  S"  t 

DISPLACEMENT  RELATIVE  TO  WANE-CAVlTN 
CROSS  SECTION  APT  or  SEPARATION  POIN  t 


Figure  2.6  Rear  axial  view  of  projectile 
cross  section. 
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Figure  2.7  Projectile  orientation  with  respect 
to  free  surface. 
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CHATTER  3 


CALCULATIONS  AND  EXPERIMENTAL  DATA 


3.1  BACKGROUND 

Detailed  information  concerning  non-normal  impact  effects  in  scarce, 
mainly  becnune  additional  parameters  that  go  beyond  those  important  for 
normal  impact  are  involved.  In  the  latter  case,  the  important  quantities 
are  the  weight,  size,  and  geometry  of  the  projectile;  the  impact  veloc¬ 
ity;  the  final  depth;  and  the  axial  deceleration  record.  For  non-normal 
teats,  however,  the  additional  important  parameters  are  the  moment  of 
inertia  and  CG  of  the  projectile,  the  lateral  components  of  the  impact 
velocity,  the  2-  or  3-D  trajectory  (including  projectile  orientation), 
and  the  lateral  components  in  the  deceleration  record. 

Usually,  the  only  way  to  obtnin  all  this  information  is  by  con¬ 
ducting  a  reverse  ballistic  test  (RBT).  Still,  if  the  projectile 
parameters  are  known,  much  insight  can  be  gained  from  knowing  only  the 
impact  conditions  and  the  final  position,  particularly  with  regard  to 
stability. 

In  this  chapter,  comparisons  will  be  made  between  calculations'*'  and 
test  results  for  RRT's  and  conventional  penetration  tests.  The  main  pur¬ 
pose  will  be  to  investigate  the  general  credibility  of  the  analysis  and 
the  parameters  used  therein.  Although  the  data  are  too  few  to  achieve 
the  degree  of  verification  that  can  be  obtained  for  normal  impact,  there 
are  still  enough  benchmarks  available  to  check  for  unreasonable  predic¬ 
tions.  With  thi3  done,  some  credence  can  be  given  to  the  parameter 
study  in  Chapter  It. 

3.2  REVERSE  BALLISTIC  TESTS  IN  SANDSTONE 

In  1077,  Sand  la  Laboratories:  conducted  four  PNA-sponsored  HRT's  In 
sandstone.  A  complete  description  of  the  tests  In  given  in  Reference  10, 


The  PENC02D  computer  code,  which  solves  the  penetration  problem  in  two 
dimensions,  using  the  analysis  presented  in  Chapter  2,  was  used  for 
all  calculations  reported  herein. 


and  the  processed  data  are  presented  in  References  11  and  12.  Two  of 
the  tests  (Nos.  2  and  3)  were  done  at  a  3-degrec  attack  angle,  and  one 
(No.  U)  at  20  degrees  obliquity  (no  angle  of  attack).  The  impact  veloc¬ 
ity  was  approximately  1500  ft/o  in  each  case.  Nominal  target  properties 


p  -  130  lb/ft3 
Y  -  3 1*00  psi 

Rock  quality  designation  (RQD)  ■  100 

A  scale  drawing  of  the  projectile  is  shown  In  Figure  3.1;  the  projectile 

parameters  are* 


hength  (L)  *  18.1*5  in 
Weight  (W)  -  <).1i8  lb 
Diameter  (D)  *  1.9  in 

Caliber  radius  head  (CRH)  ■  6.00  (beveled  tip) 

T  -  221.0  lb- in* 


For  each  case,  the  load  distributions  on  the  projectile,  calculated 

using  the  FKNC02D  computer  code  (with  output  modifications),  have  been 

used  as  input  to  a  dynamic  structural-response  code  (WHAM!')  by 

q 

T.  Uelytschko.  Predictions  were  generated  for  accelerometer  and  strain- 
gage  outputs,  samples  of  which  are  compared  with  test,  results  in 
Figures  3.2  through  3.7.  In  the  WHAM!'  on.loulat  ions ,  the  accelerometers 
indicated  in  Figures  3.2  and  3.3  were  located  on  the  axis,  l6.3'>  Inches 
from  the  nose  tip.  The  strain  gage  was  on  the  outside 

a— — - - 

The  units  shown  were  chosen  for  convenience  of  expression.  In  cal¬ 
culations  they  must  be  converted  to  a  compatible  system  (e.g,,  slugs, 

^  feet,  seconds). 

Personal  communication,  received  March  l 070,  from  T.  belytsohko,  Inc., 
Chicago,  Ill.,  to  R.  Hernard,  (Jeomeohanles  Division,  structures 
laboratory. 
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bottom  (♦“»).  9*6  inches  from  the  nose  tip. 

For  the  3-degree  attack  angle  calculation,  there  ia  clearly  a 
phase  difference  between  the  calculated  and  measured  lateral  accelera¬ 
tions  (Figure  3.2),  the  cause  of  which  is  not  yet  understood.  Also,  the 
initial  positive  peaks  in  the  calculated  results  are  somewhat  low,  but 
the  later  negative  peaks  are  in  fair  agreement  with  the  test  data. 

Better  agreement  exists  for  the  axial  decelerations  (Figure  3.3).  Fair 
agreement  is  obtained  for  the  strain  gage  output  (Figure  3.*0,  although 
the  initial  negative  strains  are  somewhat  underpredicted. 

For  the  20-degree  obliquity  calculation,  the  calculated  results  are 
in  much  better  agreement  with  the  test  data  (Figures  3.5  through  3.7), 
although  the  phase-difference  problem  is  still  apparent  in  the  lateral 
acceleration  (Figure  3.5)  data  and  the  strain  gage  output  (Figure  3.7). 

In  using  the  PENC02D  computer  code  to  generate  the  input  for 

Belytschko's  WHAMS  calculations,  the  free-surface  parameter  was  set  at 

k  =  3  ,  and  the  wake  separation  angle  was  set  at  $  .  ®  0  .  These 

min 

values  were  chosen  after  examining  the  effect  of  varying  k  and 
in  WES  predictions  for  total  lateral  force  and  pitch  moment  (Figures 
3.8  through  3.13).  The  value  k  *  8  i3  reasonable  in  light  of  pene¬ 
tration-crater  measurements  in  sandstone  (Reference  °),  and  the  A  ,  =  0 

min 

value  is  reasonable,  since  the  wake  separation  angle  should  be  small  in 
hard  targets. 

In  order  to  examine  the  effects  of  attack  angle  and  obliquity  on 
lateral  load  and  pitch  moment  in  the  WE S  theory,  a  scries  of  calcvilat ions 
has  been  made  using  the  3-  and  20-degree  RBT's  an  baseline  cases.  The 
results  are  3hovn  in  detail  in  Figures  3 . 1}*  through  3.17  and  summarised 
in  Figures  3.18  and  3.19.  A  similar  series  of  calculations  vins  made, 
varying  the  target  strength  for  the  3-degrce  RUT.  The  results,  shown 
in  Figures  3-20  and  3.21,  indicate  that  the  initial  peaks  in  the  lateral 
force  (at  0.12  ms)  are  much  more  dependent  on  strength  than  are  the 
later  peaks  (at  0.7  to  0.8  ms).  The  reason  for  the  similarity  in  the 

K 

Variation  of  k  has  no  discernible  effect  on  the  3-degree  RUT 
calculation. 


later  peaks  can  be  seen  by  examination  of  the  lateral  force  distribution 
(per  unit  length)  at  0.7  ma  (Figure  3.22).  The  amplitude  of  the  force 
distribution  1b  a  strong  function  of  the  strength,  but  the  positive  and 
negative  contributions  are  such  that  there  is  little  difference  in  the 
net  lateral  force  at  0.7  ms. 

3.3  FULL-SCALE  PENETRATION  TESTS  IN  SOIL 

Sandia  Laboratories  has  conducted  at  least  four  full-scale  3oil 
penetration  tests  that  are  suitable  as  benchmarks  for  the  WES  2-D  pene¬ 
tration  theory.  TVo  of  these  tests,  Nos.  R800915  and  R800916  (Refer¬ 
ence  13),  were  conducted  in  the  Pedro  dry  lake  site  on  Tonopah  Test 
Range  (TTR).  The  other  tests,  Nob.  Rb5b025-22  and  Rb5b025-23  (Refer¬ 
ence  lU ) ,  took  place  in  the  Antelope  dry  lake  site  on  TTR.  Pertinent 
information  for  all  four  tests  is  summarized  in  Table  3.1,  and  Young'3 
S-number  profile  (Reference  lb)  for  the  Antelope  site  is  given  in 
Table  3.2. 

Using  the  information  given  in  Tables  3-1  and  3.2  as  input, 
calculations  have  been  made  for  the  penetrator  trajectories.  The  re¬ 
sults  are  compared  vith  the  Sandia  test  data  in  Figures  3.23  through 
3.27. 

Calculations  for  teats  R800915  and  R800916,  presented  respectively 
in  Figures  3.23  and  3.2b,  show  the  effects  of  varying  the  wake  separa¬ 
tion  angle.  Here  the  experimental  results  are  inconsistent.  The  tests 
were  Identical  except  for  a  20  percent  difference  in  impact  velocity, 

yet  there  was  a  500  percent  difference  in  lateral  (horizontal)  displace- 

5 

ment.  One  of  the  reasons  for  the  discrepancy  may  have  been  a  differ¬ 

ence  in  the  amount  of  wake  separation  In  the  two  tests.  Separation 
alone,  however,  in  not  enough  to  account  for  the  difference,  since  even 
the  calculation  for  $  ^  *  0  predicts  3  feet  of  lateral  displacement 
for  test  R8OOOI5.  Tn  these  calculations  the  free-surface  parameter  was 


The  average  vertical  penetration  resistance  apparently  was  the  same, 
because  the  depths  can  be  calculated  with  same  2-number . 


« 
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set  at  k  »  0  ,  due  to  its  lack  of  influence  on  near-vertical  impact 
problems. 

Test  R*t5**025-23  is  a  good  example  for  assessing  the  effects  of  the 
free*-surface  parameter  independent  of  the  vake  separation  angle  effects. 
Here  the  initial  obliquity  is  large,  and  k  and  $min  both  influence 
the  calculated  trajectory  (Figures  3.25  and  3.26).  For  this  test  the 
best  simulation  of  the  observed  trajectory  is  obtained  for  k  *  *♦  and 
<t  .  *  0  Nevertheless,  it  is  clear  from  the  calculations  that  the 

vake  separation  angle  has  a  much  stronger  influence  than  the  free- 
surface  parameter,  even  for  large  obliquity. 

Test  R**5*»025-22  was  conducted  with  a  shorter  projectile  (L/D  ■  8) 
than  the  other  three  tests  (L/D  *  10).  In  this  case,  the  penetrator 
impacted  at  30  degrees  and  rotated  80  degrees  before  coming  to  rest  13-5 
feet  below  the  target  surface.  With  the  free-surface  parameter  set  at 
k  *  U  ,  calculations  made  for  different  values  of  the  vake  separation 
angle  (Figure  3.27)  indicate  that  a  value  of  $  ^  a  *+  degrees  best 
reproduces  the  observed  final  position  of  the  penetrator. 

In  other  (normal  impact)  tests  in  the  Antelope  site,  Sandia  has 
found  L/D  =  10  to  be  significantly  more  stable  than  L/D  *  8  (Refer¬ 
ence  1*0.  For  comparison  of  the  two  projectiles  in  this  study,  calcula¬ 
tions  were  made  for  the  shorter  projectile  (from  test  R*<5*<025-22 )  using 
the  S-number  and  impact  conditions  of  test  R800916,  setting  k  =  0  ,  and 

varying  $  .  .  The  results  (Figures  3.28-3.31)  indicate  that  the 
min 

stability  of  either  projectile  is  strongly  dependent  on  the  value 
selected  for  the  vake  separation  angle. 


To  simulate  the  test  conditions  (Table  3.1),  the  calculation  was 
_  started  with  the  nose  tip  at  Z  =  8  inches  . 

8=2  degrees,  a  <  0,5  degrees. 


Table  3.2.  S-number  profile  for  TTR 
Antelope  Dry  Lake  site.a 
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0 
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7 

; 

8o 

85 

3 

> 

85 
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STATION 

AXIAL 
DISTANCE 
FROM 
NOSE  TIP 
INCHES 

SIGNIFICANCE 

a 

3.80 

NOSE  BASE 

CG 

8.71 

CENTER  OF  GRAVITY 

F 

8.20 

FORWARD  ACCELEROMETER 

R 

17.30 

REAR  ACCELEROMETER 

S 

9.60 

STRAIN  GAGE 

NOTE:  FLARE  BEGINS  AT  S.64IN  AND 
ENOS  AT  13.04  IN 
TOTAL  LENGTH  L  3  |«.  IS  IN 

CG  B 


~rf 

R  3  F 


I _ I - 1 - 1 — 

20  IS  10  S 

DISTANCE  FROM  NOSE  TIP,  INCHES 


Figure  3.1  Scale  drawing  of  RBT  penetrator 
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Figure  3.2  Comparison  of  calculated  and 
measured  lateral  accelera¬ 
tions  for  Sandia/DNA  RBT, 
3-degree  attack  angle. 


Figure  3.3  Comparison  of  calculated  and 
measured  axial  decelerations 
for  Sandia/DNA  RBT,  3-degree 
attack  angle. 
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Figure  3.6  Comparison  of  calculated  and  mea¬ 
sured  axial  decelerations  for 
Sundia/PNA  RBT,  20-degree 
obliquity. 


01  0  4  0*  0  I  10 

T  t  M •  ,  MU 

Figure  3.7  Comparison  of  calculated  and  mea¬ 
sured  strains  for  Band i a /ON A  KPT, 
20-degree  obliquity. 
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Figure  3.22  Calculated  lateral 
force  distribution 
for  different  values 
of  target  strength, 
Sandia/DHA  RBT,  3- 
degree  attach  angle. 


Figure  3.23  Effect  of  wake 
separation 
angle  $min 
on  calculated 
trajectories , 
Sandia  Test 
R800915. 
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Figure  3.27  Effect  of  wake  separation 

angle  *min  on  calculated 
trajectories,  Sandia  Test 
Rl*  5^025-22. 


MO«tt*NT«V  MltTtON,  *▼ 


V«*TlCAi-  .  0T 


Figure  3.30 


Comparison  of  calculated  tra¬ 
jectories  for  Sandia  pene- 
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Figure  3.31  Comparison  of  calculated  tra¬ 
jectories  for  Sandia  pene- 
trators,  4m^n  *  8  degrees. 
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CHAPTER  4 


PARAMETER  STUDY  FOR  SOIL  PENETRATION 

4.1  BACKGROUND 

The  PENC02D  computer  code  predicts  trajectories  that  seem  reasonable 
in  light  of  existing  non-normal  soil  penetration  data  (Chapter  3)* 
although  inconsistencies  in  the  test  results  leave  unresolved  the  ques¬ 
tion  of  vake  separation  and  its  quantitative  effect  on  stability.  Un¬ 
certainties  notwithstanding,  enough  benchmarks  do  exist  so  that  some 
credence  can  be  given  to  a  study  of  the  effects  of  varying  the  projectile 
parameters,  the  impact  conditions,  and  the  target  penetrability. 

A  baseline  projectile,  target,  and  set  of  impact  conditions  are 
specified  in  this  chapter,  and  calculations  are  presented  in  which  indi¬ 
vidual  parameters  have  been  varied  one  at  a  time.  In  some  cases  the 
variation  of  a  given  quantity  independent  of  other  quantities  may  be  un¬ 
realistic  in  a  practical  sense  (e.g.,  changing  the  projectile  length 
without  affecting  the  moment  of  inertia).  However,  the  objective  in 
such  instances  is  to  show  the  influence  of  a  particular  parameter  in  the 
calculation  itself.  In  a  practical  design-parameter  study  (Section  4.6), 
coupled  parameters  (e.g.,  weight,  length,  and  mass  moment  of  inertia) 
have  to  be  varied  together. 

4.2  BASELINE  CONDITIONS 

The  projectile  chosen  for  this  study  is  similar  to  the  one  used  in 
Sandia  Test  R454025-23  (L/D  »  8)  and  represents  a  marginally  stable 
design,  according  to  Sandia's  experience.  The  pertinent  projectile 
quantities  are 

W  a  300  lb 
D  *  6  in 
L  *  48  in 
CRH  »  6.o 

CG  location  c^o  *  26.4  in  from  nose  tip 


.  U  2 
I  *  U  x  10  lb-in 

Aftbody  taper  angle  $  ^  *  0 

The  baseline  impact  conditions  are 

V  *  1500  ft/sec 
o 

a  U  degrees 
9q  *  4  degrees 

The  baseline  target  properties  are 

S  *  8 
k  a  U 

4 

min  ~  2  degrees 

All  of  the  above  parameters  will  be  varied  except  k  and  $  .  ,  whose 

effects  have  already  been  demonstrated  in  Chapter  3.  Aside  from  the 
individual  quantity  or  quantities  being  varied,  the  input  will  be  the 
same  as  the  baseline  conditions.1 

Figures  U.l,  U.2,  and  U.3  show  the  trajectory,  external  forces,  and 
pitch  moment,  respectively,  calculated  for  the  baseline  conditions. 

Since  the  baseline  calculation  uses  a  constant  wake  separation  angle 
^min  "  ^  degrees),  the  penetration  path  is  curved,  rather  than  straight, 
and  there  is  some  residual  oscillation  in  the  pitch  moment  as  aftbody 
contact  alternately  increases  and  decreases. 

U.3  VARIATION  OF  TARGET  PENETRABILITY 
AND  IMPACT  CONDITIONS 

Figure  U.U  indicates  the  degree  to  which  the  target  S-number  in¬ 
fluences  the  trajectory.  The  main  effect  seems  to  be  an  increase  in 
the  total  path  length  with  increasing  S  .  Although  the  horizontal  dis¬ 
placement  increases  with  S  ,  there  seems  to  be  little  variation  in  the 

1  The  baseline  calculation  will  appear  in  all  comparisons  of  trajectories 
except  Figure  U.7,  and  will  be  represented  by  a  solid  curve  (e.g., 
Figure  U.l) . 
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overall  path  shape.  The  same  thing  can  be  said  for  the  impact  velocity 
(Figure  U . 5 ) . 

The  initial  attack  angle  has  a  mild  effect  on  the  trajectory  (Fig¬ 
ure  4.6),  aside  from  a  marked  increase  in  path  curvature  and  horizontal 
displacement  between  0  and  2  degrees.  The  Q-degree  case  goes  straight 
in  because  there  is  no  angle  of  attack  (or  obliquity)  to  initiate  rota¬ 
tion.  Large  obliquities  (oq  *  0,  >  60  degrees)  can  produce  ricochet 

(Figure  4.7),  although  lesser  obliquities  seem  to  have  a  gentle  effect 
on  path  curvature.2 

U . U  INDEPENDENT  VARIATION  OF  PROJECTILE 

PARAMETERS  FOR  STRAIGHT  AFTBODY 

Figures  4.8,  4.9,  and  4.10  show  the  effects  of  varying  the  weight, 
diameter,  and  nose  shape,  respectively.  Treated  as  independent  param¬ 
eters,  none  of  these  variables  has  much  effect  other  than  changing  the 
path  length.  The  horizontal  displacement  increases  with  path  length, 
but  the  shape  of  the  trajectory  changes  only  slightly. 

It  might  be  argued  that  the  mass  moment  of  inertia  alone  should 
have  a  big  influence  on  the  trajectory.  This  is  not  the  case,  however, 
as  can  be  seen  in  Figure  4.11.  Although  the  mass  moment  of  inertia  ha3 
a  moderately  stabilizing  effect,  it  is  overshadowed  by  aftbody  separation, 
which  allows  the  projectile  to  rotate  with  little  opposition  during  most 
of  the  event  (see  Figure  4.3). 

Figure  4.12  shows  the  consequences  of  treating  the  projectile  length 
as  an  Independent  parameter,  keeping  the  geometric  position  of  the  CG 
constant  ( CC(./L  =  0.55).  A  reduction  in  projectile  length  has  a  stabi¬ 
lizing  effect,  due  to  reductions  in  the  moment  arm  and  the  total  (lateral) 
loaded  area  on  either  side  of  the  CG.  This  is  the  least  realistic  of 
all  the  foregoing  trajectory  comparisons  because  a  physical  change  in 
length  changes  all  other  parameters  (Section  4.6). 

The  most  important  single  parameter  is  the  geometric  CG  location 


The  initial  obliquity  required  for  ricochet  would  be  smaller  if  a 

larger  value  has  been  chosen  for  $  , 

min 


48 


X 


C  _/L  ,  as  is  indicated  in  Figure  U.13.  If  the  CG  is  far  enough  forward 
(5„„/L  ■  0.1*5),  the  negative  moment  (after  CG  entry  into  the  target) 
outweighs  the  initial  positive  moment,  and  a  rotation  reversal  occurs. 

If  the  CG  is  far  enough  aft  (cc(,/L  3  0.65),  the  initial  positive  moment 
is  so  great  that  the  projectile  ends  up  going  sideways. 

U.5  VARIATION  OF  AFTBODY  FLARE 

Figure  U.lU  shows  trajectories  calculated  for  different  values  of 

the  aftbody  flare  angle  .  In  each  calculation  the  aftbody  merges 

with  the  ogive  nose  section  at  the  point  of  tangency,  where  n  =  • 

An  illustration  for  $  3  3  degrees  is  shown  in  Figure  l*. 15. 

a  rt 

If  *  „  =  a  ,  there  is  no  loss  of  contact  with  the  target,  and 
aft  min 

the  penetrator  goes  almost  straight  in.  On  the  other  hand,  if 

0  <$_.<$  .  ,  the  trajectory  is  straighter  than  the  baseline  case 

art  min 

(straight  aftbody,  $  *  0)  but  it  is  still  noticeably  curved.  Flare 

ax  r 

angles  greater  than  $  .  reduce  the  final  depth,  due  to  increased  base 

min 

diameter  (decreased  W/A). 

From  this  comparison,  the  optimum  flare  angle  appears  to  be  a  value 
less  than  the  wake  separation  angle,  producing  a  slightly  curved  trajec¬ 
tory  but  leaving  the  depth  essentially  unaffected.  The  flare  must  be 
installed  without  shifting  the  CG  toward  the  tail. 

U.6  VARIATION  OF  COUPLED  PROJECTILE  PARAMETERS 

The  projectile  parameters  that  were  varied  independently  in  Sec¬ 
tion  4.U  are  coupled  to  one  another  in  practice.  As  far  as  stability  is 
concerned,  the  most  important  coupling  occurs  among  weight,  length,  CG 
location,  and  mass  moment  of  inertia.  A  serious  parameter  study  should 
consider  the  simultaneous  variation  of  these  quantities. 

A  simple  but  fairly  realistic  way  of  accounting  for  the  change  in 

W  ,  £  ,  and  I  with  L  is  to  add  or  subtract  cylindrical  aftbody 

CG 

sections  of  constant  linear  density  from  the  baseline  projectile,  while 


3 


Constant  weight/unit  length  =  300  lb/H8  in  =  6.25  lb/in. 


I 


2 

assuming  I  *  ML  .  Figure  4.16  shove  trajectories  for  projectiles 

whose  properties  vere  varied  in  this  way. 

According  to  the  calculations,  if  the  projectile  is  short  enougn 

(L/D  <  5),  the  CG  will  lie  in  an  unstable  position,  ultimately  producing 

U 

sidevays  motion.  Otherwise,  the  effect  of  the  changing  parameters  is 
little  more  than  an  increase  in  path  length  with  projectile  length. 


See  also  Figure  U.13, 
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Figure  U.l  Calculated  trajectory  for 
baneline  conditions. 
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Figure  I*. 2  Calculated  external  forces  for  baseline 
conditions . 
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Figure  U.3  Calculated  pitch  moment  for  bnBelitie 
conditions. 
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Figure  >».!i  Effect  of  soil  penetrability 
on  calculated  trajectories. 
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Figure  U.15 


Scale  drawing  of  projectile  with  3- 
degree  aftbody  flare 


Figure  *».l6 


Combined  effects  of  weight, 
length,  CG  location,  and 
mass  moment  of  inertia  on 
calculated  trajectories. 


CHAPTER  5 


CONCLUSIONS 

An  empirical  theory  has  been  developed  for  analyzing  impact  and 
penetration  in  soil  and  rock.  Calculations  with  this  theory  reproduce 
penetration  data  for  normal  impact  at  least  as  well  as  other  empirical 
analyses  (e.g.,  Young's  equation).  The  extension  from  one-  to  two- 
dimensional  motion  was  obtained  by  extrapolation,  with  simple  models 
postulated  for  the  frce-surface  and  wake  separation  effects.  The  entire 
non-normal  analysis  was  then  incorporated  into  the  PENC02D  computer  code. 
Benchmarks  for  the  theory  consist  of  strain  and  acceleration  data  from 
RBT's  in  sandstone,  and  trajectory  data  from  conventional  penetration 
tests  into  soil. 

Projectile  acceleration  and  strain  predictions,^"  generated  from  WES- 
calculated  external  loads,  agree  fairly  well  with  the  RBT  data.  There  is, 
however,  a  tendency  to  underpredict  the  lateral  loads  for  yawed  impact 
somewhat.  There  is  also  an  unexplained  phase  difference  between  the  cal¬ 
culated  and  observed  lateral  accelerations  for  both  yawed  and  non-yawed 
oblique  impact.  On  the  other  hand,  the  calculated  variations  of  maximum 
lateral  force  and  pitch  moment  with  initial  obliquity  and  attack  angle 
are  quite  reasonable  (essentially  linear).  In  any  case,  the  recommended 

values  for  the  wake  separation  angle  and  free-surface  parameter,  respec- 

o 

tively,  are  *  0  and  k  »  8  for  rock  and  rocklike  materials. 

Comparisons  of  calculated  and  observed  penetrator  trajectories  in 
soil  are  inconclusive,  since  the  test  data  are  not  consistent.  Concern¬ 
ing  terradynamic  stability,  however,  the  most  important  target  parameter 

3 

is  the  wake  separation  angle,  whereas  the  most  important  projectile 


WHAMS  code  (dynamic  structural  response)  calculations  made  by 
2  T.  Belytschko,  with  PENC02D  calculations  as  input. 

^  These  values  gave  reasonable  external  load  predictions  for  the  RBT's. 

At  this  point,  it  is  impossible  to  specify  the  most  likely  values  of 
the  free-surface  parameter  and  the  vake  separation  angle,  A  reasonable 
guess  for  the  respective  ranges  of  values,  however,  is  0  k  <  3 
and  0  <_  4  <8  degrees  . 
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parameter  is  the  CG  location.  Aftbody  flare  also  has  a  beneficial 
effect,  if  the  CG  is  not  shifted  toward  the  tail  and  the  base  diameter 
is  not  increased  enough  to  significantly  reduce  the  final  depth. 

One  phenomenon  not  predicted  by  the  present  theory  is  the  marked 
improvement  in  projectile  stability  achieved  by  going  from  L/D  ■  8  to 
L/D  ■  10  ,  which  was  observed  in  tests  conducted  by  Sandia  Laboratories 
(Reference  lU).  WES  calculations  indicate  a  gradual  improvement  in 
stability  with  projectile  length  (see  Figure  U.l6).  If  L/D  *  8  does 
represent  a  critical  value  for  projectile  stability,  this  vould  indicate 
that  some  basic  mechanism  is  missing  in  the  theory,  either  in  the  wake 
separation  analysis  or  in  the  stress  equations. 

One  mechanism  conspicuously  absent  in  the  present  wake  separation 

analysis  is  the  resistance  of  the  target  to  the  (cylindrical)  expansion 

of  the  wake  cavity,  which  would  lead  ultimately  to  the  rebound/collapse 

of  the  cavity  on  the  aftbody.  The  probable  variation  of  4  .  with 

sun 

velocity,  depth,  nose  shape,  and  target  penetrability  (S-number)  is  also 
neglected.  Finally,  the  stress  equations  may  be  underpredicting  the 
effects  of  lateral  motion.  Future  work  should  consider  these  problems 
and  any  others  that  might  bear  on  the  stability  question. 


fvjav 
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The  latter  conclusion  was  also  drawn  from  the  preliminary  analysis 
(Reference  8). 
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APPENDIX  A 


PENETRATION  THEORY  FOR  ROCK  AND  CONCRETE 

A.l  EMPIRICAL  THEORY  FOR  NORMAL  IMPACT 

Early  penetration  analyses  at  WES  (References  15  through  IT)  relied 
heavily  on  the  Cavity  Expansion  Theory  (CET),  whereby  the  coefficients 
a  and  b  and  the  function  f(V)  in  the  resisting  stress 

a  *  a  +  b  f(V)  (A.l) 

could  be  obtained  explicitly  in  terms  of  the  density,  strength,  and 
elastic  moduli  of  the  target.  Examination  of  penetration  data  from 
several  sources,  however,  has  failed  to  show  a  consistent  correlation 
between  penetrability  and  standard  engineering  properties,  especially 
for  soil.  Although  the  most  consistently  accurate  predictions  have 
been  obtained  for  concrete,  additional  factors  such  as  aggregate  3ize 
may  overshadow  elastic  effects. 

In  the  long  run,  a  " fundament al"  theory  like  the  CET  requires 
empirical  modification  in  order  to  be  of  general  use  as  a  predictive 
tool.  An  empirical  approach  from  the  outset  seems  better,  thus  keeping 
the  functional  relations  among  the  parameters  as  simple  as  possible. 

For  concrete  and  rock,  a  stress  equation  of  the  form 

a  =  a  +  bV  (A. 2) 

does  seem  adequate  for  correlating  penetration  data  obtained  with  a 
given  projectile/target  combination.  In  general,  however,  the  coeffi¬ 
cients  a  and  b  will  vary  with  the  characteristics  of  the  projectile 
and  the  target.  Based  on  penetration  tests  in  concrete  (Reference  18) 
and  intact  rock  (Reference  19),  the  following  empirical  formula  is  pro¬ 
posed  for  the  stress  distribution  on  the  projectile1 

1  All  dimensional  quantities  must  be  expressed  in  compatible  units; 
e.g.,  slugs ,  feet,  and  seconds. 


6  5 


(A.  3) 


o  -  1.!>8P  (>*Y  ♦  IvSfrt)  Ain  n 


Relating  Equation  A.  3  to  A .P,  it  in  clear  that  a  ■  6.3P8  Y  /nin  n 
and  thnt  b  ■  k.J^6  /pY  nin  n  .  When  Equation  A. 3  in  nubntituted  into 
Equations*  P.lt  and  P.8,  the  equation  of  motion  (Equation  P.10,  0  ■  0) 

for  n  fully  embedded  projectile  becomes) 

MZ  •  (  >tY  ♦  3V^p7)  (A. I*) 

"re 

where  D  in  the  maximum  diameter  and  N^c  in  the  none  performance 
coefficient  in  rock  and  concrete,  given  for  ogiven  by 
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rc 


O.863 


ItOHIl  -  1 
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and  for  conen  by 
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0.808 
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(A. 6) 


with  n,  reprenentlng  the  cone  half-angle  ami  (’811  reprenent.ing  the 
rndiun  of  curvature  of  the  tangent  ogive  in  cnllbern  (multipten  of  P) . 

Pi  nregardlng  tlie  changing  cronn  nee  1. 1  on  during  the  none-embedment 
procenn,  the  final-depth  no  .tut.  ion  to  Equation  A.  It  in 
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(A. 7) 


wltere  l:t  the  impact,  velocity. 

In  mid  it  loti  to  the  parameters'  given  in  Equation  A.  7,  there  in  at 
leant  one  more  parameter  that  ir.fiuencess  penet  rabl  1  i  t  y .  In  “ook  tlii:i 
quantity  in  the  Kook  Quality  Pen  I  gnat,  ion  (Ks,'P),  introduced  hy  Peere 
(Reference  P0).  The  RiiP,  an  index  for  the  degree  of  in  nlt.u  fracturing 


t>(> 


in  a  given  site,  is  obtained  using  a  modified  core-logging  procedure: 

All  solid  pieces  of  core  that  are  4  inches  long  or 
longer  are  added  up,  and  this  length  is  called  the  mod¬ 
ified  core  recovery.  The  modified  core  recovery  is  di¬ 
vided  by  the  total  length  of  core  run,  and  the  quotient 
multiplied  by  100  percent  is  the  value  of  the  RQD. 

In  concrete,  the  maximum  size  of  coarse  aggregate  D  apparently  has  a 

ag 

mild  effect  on  penetrability  (Reference  21),  via  the  ratio  D  /D  . 

ag 

Letting  Y  represent  the  Intact  unconfined  strength  for  rock  and 
concrete,  the  effective  strengths  for  penetration  can  be  calculated  from 


,r  ■ ,  •  ©r  u.8) 

/3D  \0-2 

h  ■ J  ’  (,-d“7  (a-9) 


Thus,  to  calculate  penetration  in  less-than-intact  rock  (RQD  <  100),  the 

quantity  Yf  should  replace  Y  in  Equations  A. 3,  A.U,  A. 7,  and  2.12. 

Likevise,  when  the  maximum  aggregate  size  is  specified  for  a  given 

2 

concrete,  Yc  should  replace  Y  in  the  same  equations. 

Figure  A.l  shows  a  comparison  of  Equation  A. 7  with  nondimension- 

alized  penetration  data  for  rock  and  concrete.  The  concrete  data  were 

obtained  for  a  case  in  which  D/D  =  3  ;  i.e.,  Y  =  Y  .  Primarily, 

ag  c 

the  figure  shows  the  degree  to  which  Equation  A. 7  is  able  to  collapse 

the  data  to  a  single  curve.  There  is  some  residual  scatter,  but  thi3 

correlation  represents  an  improvement  over  previous  attempts  (References 

...  3.  .  '  ••  . . 

U  and  22).  Table  A.l  presents  the  data  in  dimensional  form,  whereas 

Table  A. 2  gives  the  target  properties  and  projectile  parameters. 

A. 2  EXTRAPOLATION  TO  NON-NORMAL  IMPACT 

Equation  A. 3  was  obtained  by  starting  with  Equation  A. 2  and 


In  structural  concrete,  Dag  usually  varies  from  0.75  to  1.5  inches, 
o  The  RQD  is  always  assumed  to  be  100  for  concrete. 

The  improvement  comes  from  the  addition  of  a  nose-shape  effect  and  the 
reduction  of  the  RQD  effect  initially  used  in  References  and  22. 
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adjusting  coefficients  until  a  good  data  fit  was  achieved  for  combined 
rock  and  concrete  penetration  data  (Figure  A.l)  taken  from  References 
18,  19,  and  23.^  Previous  experience  with  these  materials  had  shown 
that  a  and  b  should  be  proportional  to  Y  and  ,  respectively; 

the  /sin  ^-proportionality  (nose-shape  effect)  waB  drawn  from  Young's 
observations  of  soil  penetration  and  was  presumed  to  hold  true  for 
harder  targets.  With  the  addition  of  the  RQD  effect  on  the  apparent 
strength,  the  correlation  shown  in  Figure  A.l  was  achieved. 

The  form  of  Equation  A. 3  is  such  that  it  can  readily  be  extrapo¬ 
lated  to  the  case  of  nonaxial  motion.  Recognizing  that  for  axial 
motion 

V 

sin  n  ■  y2-  (A. 10) 


where  Vn  =  component  of  V  normal  to  projectile  surface,  it  i3  reason¬ 
able  to  replace  V  and  V  ,  respectively,  by  the  local  velocity  v 

n 

and  its  normal  component  vn  .  Thus,  for  2-  or  3-dimensional  motion, 
the  expression  for  the  compressive  normal  stress  on  the  projectile 
surface  is 


a  * 


1.582  (UY  +  3v/pY 


(A. 11) 


where 

v  *  absolute  local  velocity  relative  to  fixed  target 
v^  »  component  of.  v  normal-  te  projectile'  surface 

In  order  for  a  stress  to  act  on  a  given  projectile  surface  element, 
there  must  be  contact  with  the  target  (Section  2.k)  and  net  motion  into 
the  target  (v  >  0).  Otherwise,  it  is  presumed  that  a  =»  0  . 


c  Additional  material  property  information  is  given  in  References  2 h— 26 . 
3  For  purely  axial  motion,  V  and  Vn  are  identical  with  v  and  vn  , 
respectively.  In  extrapolating  to  a  local  definition  of  stress  for 
nonaxial  motion,  the  most  direct  route  is  to  use  the  local  velocity. 
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Table  A.l.  Penetration  data  for  rock  and  concrete 


Test  Source _ 

Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 
Canfield  and  Clator 

Sandia/DNA 

Sandia/DNA 

Sandia/DNA 

Sandia/DNA 

Sandia/DNA 

Sandia/DNA 

Sandia/DNA 

Sandia  No.  120-77 

Sandia  No.  120-112 

Sandia  No.  120-103 

Sandia  No.  120-106 


Target  Material 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Concrete 

Welded  Tuff 
Welded  Tuff 
Welded  Tuff 
Welded  Tuff 
Welded  Tuff 

Sandstone 

Sandstone 

Welded  Agglomerate 
Sandstone 
Sandstone 
Granite 


Impact 

Velocity 

fps 

Penetrat: 

Depth 

in 

1005 

8.0 

1025 

9.0 

1250 

10.0 

11*85 

H*. 5 

1775 

16.5 

1975 

23.5 

2020 

19.5 

2325 

26.0 

2350 

21+  .0 

21*30 

27.5 

2535 

29.0 

2655 

29.5 

1220 

87 

1350 

102 

1560 

ll*2 

161*0 

132 

1650- 

132 

ll*55 

ll*0 

1505 

ll*6 

1065 

156 

821* 

122 

880 

120a 

860 

150 

In  this  test  the  sandstone  was  covered  by  30  inches  of  soil,  so  the 
toted  depth  of  penetration  is  150  inches.  However,  the  penetration 
resistance  of  the  soil  is  negligible  compared  with  that  of  the  rock. 


( 
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Table  A.  2.  Projectile  and  target  parameters. 
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GRANITE 


SOURCE  OF  DATA 

CANFIELD  &  CLATOR 
SAN  01  A/ON  A  TESTS 
SANOIA/ONA  TESTS 
SAN 01 A  TEST  NO.  120-77 
SAN 01 A  TEST  NO.  120-  1 12 
SANOIA  TEST  NO.  120-  103 
SAN 01 A  TEST  NO.  120-  106 


Figure  A.l  Nondimen sionali zed  penetration  data  for  rock 
and  concrete. 
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PENETRATION  THEORY  FOR  SOIL 

B.l  EMPIRICAL  THEORY  FOR  NORMAL  IMPACT 

To  date.  Young's  equation  (References  1-3)  has  been  the  most  con¬ 
sistently  successful  formula  for  calculating  penetration  depths  in  soil. 
Rather  than  using  standard  engineering  properties  to  describe  the 
target.  Young  lumps  the  penetrability  into  a  single  empirical  variable, 
represented  by  S  .  In  U,  S.  customary  units Young's  equation  is 

0.0031  S  Ng  ft]k  (vq  -  IOO)  ,  VQ  >_  200  ft/s 

___  /  2V  2 

0.53  3  »3  m  to  \l  *  3^07000 

where 

Z  =  final  depth,  ft 
max 

W  =  projectile  weight,  lb 

7T  2 

A*  jD  =  maximum  cross-sectional  area,  sq  in 
Vq  =  impact  velocity,  ft/s 

Ng  =  nose  performance  coefficient  in  soil,  dimensionless 
S  =  soil  penetrability  index,  dimensionless 
The  S-number  can  be  obtained  directly  from  previous  penetration 
data  for  a  given  target,  or  it  can  be  estimated  from  data  for  geologi¬ 
cally  similar  target's  (Table  B.l/.  '  Whatever  the  case.  Young's  equation 
represents  a  scaling  relation  between  the  projectile  characteristics 
(size,  weight,  geometry),  the  impact  velocity,  and  the  maximum  depth  of 
penetration. 

Since  Young's  equation  gives  only  the  final  depth,  it  cannot  be 
used  to  calculate  the  instantaneous  resisting  force,  only  the  average 
resisting  force.  This  is  not  a  drawback  for  normal  impact,  since  the 


Young's  equation  is  not  dimensionally  homogeneous.  All  quantities 
must  be  specified  in  the  units  given. 


j  ,  Vq  <  200  ft/s 


axial  resisting  force  is  often  a  step-pulse  (i.e.,  uniform*  except  for 
some  oscillation)  in  a  given  layer.  On  the  other  hand*  for  non-normal 
impact,  both  the  axial  and  lateral  forces  must  he  known  if  the  2-  and 
3-D  motion  is  to  he  calculated.  Since  Young's  equation  gives  no  indica¬ 
tion  of  the  lateral  forces,  a  more  detailed  analysis  is  needed. 

In  Reference  U  a  projectile  equation  of  motion  vas  developed  for 
normal  impact  and  penetration  in  soil,  based  on  data  originally  used 
by  Young.  Postulating  that  the  axial  force  should  (a)  have  the  fora  of 
a  step-pulse  and  (b)  exhibit  linear  dependence  on  velocity,  it  was 
found  that  an  equation  of  the  form 

-M  ~  »  a  +  bV  +  cZ  (B.l) 

at 

would  satisfy  "both  these  criteria*  provided  that 

b  «  /Ste  (B.2) 

The  explicit  appearance  of  the  projectile  mass  in  the  velocity  coeffi¬ 
cient  b  was  puzzling  but  absolutely  necessary  to  generate  a  step- 
pulse  resisting  force.  Moreover,  it  brought  the  mass  dependence  of  the 
maximum  depth  into  agreement  with  Young.  Physically,  however,  the  mass 
should  occur  only  on  the  left  side  of  Equation  B.l,  with  parameters  such 
as  diameter,  nose  shape,  and  S-number  appearing  in  the  coefficients  a  , 
b  ,  and  c  . 

The  mass  dependence  of  the  coefficient  b  is  now  thought  to  be  a 

3ize  effect  in  disguise.  This  is  reasonable  since  many  of  the  larger 

3  2 

projectiles  used  by  Young  had  masses  roughly  proportional  to  D  . 
Consequently,  an  alternative  to  Equation  B.2  is 


The  deep  earth  penetrators,  built  and  tested  by  Sandia  Laboratories, 
had  L/D  ratios  of  about  10  to  maintain  stability.  For  many  of  these 
projectiles  the  weight  (in  pounds)  can  be  related  to  the  maximum 
diameter  (in  inches)  by 


W  =  1.3  D3  +  30* 


7U 


2 


(B.3) 


and  the  dependence  of  a  and  c  upon  S  and  D  is  the  same  as  in 
Reference  U: 

(B.5) 

(B.6) 

which  insures  that  for  low-velocity,  shallow  penetration. 


in  accordance  with  Young's  equation.  When  the  nose  performance  coeffi¬ 
cient  Ng  is  added,  the  coefficients  a  ,  b  ,  and  c  can  be  expressed 
in  terms  of  three  "universal"  constants  u  ,  &  ,  y  : 


a 


S  N 

3 


(B.8) 

(B.9) 

(B.10) 


With  these  expressions  for  a  ,  b  ,  and  c  ,  it  follows  from  Equa¬ 
tion  B.l*  that 
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where 


1  <  n  <  2.5 


(B.12) 


and  Young's  nose  performance  coefficient  in  soil  is  approximated  by 


0.632 


U(crh)* 

i*  CRH  - 


1/U 


ogives 


(B.13) 


,  cones  (B.l4) 

Vain  nc 

The  values  of  the  constants  p  ,  0  ,  and  y  have  been  obtained. 

The  ratio  of  0/y  was  adjusted  until  step-pulse  deceleration  curves 
were  achieved.  The  values  of  u  and  y  were  then  determined  by 
fitting  Equation  B.4  to  shallow  penetration  data  (Reference  27)  for 
projectiles  with  the  same  weight  and  diameter^  in  the  Main  Lake  area  of 
Tonopah  Test  Range  (TTR),  The  values  for  the  three  constants*1  are 

u  *  1.96  *  105  lb/ft 
0  »  6.25  x  106  lb2s2/ftT 
Y  *  1.56  x  10?  lb/ft 


Figure  B.l  compares  Equation  B.4  with  normalized"5  deep-  and 


3 

£  W  ■  22U  pounds,  D  *  4. 4  inches  . 

Equations  B.l,  B.4,  B.8-B.10,  B.15,  and  B.l6  are  dimensionally 
homogeneous.  Any  set  of  dimensionally  compatible  units  may  be  used 
therein,  as  long  as  all  quantities  are  converted  to  the  same  system 
of  units. 

'  Normalized  using  Equation  B.4.  Amount  of  scatter  in  normalized  data 
indicates  degree  of  validity  as  scaling  relation. 
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shallow-penetration  data  for  TTR  Main  Lake.  The  same  data  are  given  in 
original  form  in  Reference  2,  where  Young  gives  the  S-number  profile  as 
S  ■  5*2  from  0  to  8  feet  and  S  5  2.5  from  8  to  25  feet.  Equation  B.4 
agrees  best  with  the  shallow  data  using  S  =  4.3  and  with  the  deep 
data  using  S  *  3  . 

Equation  B.4  appears  to  do  about  as  well  as  Young's  equation  in 
extrapolating  and  interpolating  soil  penetration  data,  but  as  a  formula 
for  depth  prediction,  it  offers  no  advantage  over  the  latter.  On  the 
other  hand,  the  underlying  equation  of  motion  (Equation  B.l)  contains 
some  additional  information  concerning  the  explicit  relation  between 
the  axial  resisting  force  and  the  velocity,  depth,  diameter,  and 
geometry  of  the  projectile. 

The  analysis  of  nonaxial  motion  requires  that  the  stress  distribu¬ 
tion  be  specified  on  the  projectile.  In  the  present  context,  any  stress 
distribution  is  potentially  correct  if  it  generates  Equation  B.4  with 
the  values  of  a  ,  b  ,  and  c  given  by  Equations  B.8-B.10.  The 
simplest  distribution  that  accomplishes  this,  giving  approximately  the 
right  dependence  of  Ng  on  projectile  geometry,  is  a  compressive  normal 
stress  defined  as  follows: 


/sin  n  + 


S  r 


S  ‘p 


(B.15) 


where 


Tp  *  local  cylindrical  radius  of  projectile 

H  *  local  angle  be^weep  projectile  surface  and  axis  of  symmetry 
‘  (Ffgure  2.1 J 

Z  *  vertical  depth  from  target  surface  to  point  on  projectile 
surface  (Figure  2.7) 


B.2  EXTRAPOLATION  TO  NON-NORMAL  IMPACT 

Now  that  the  stress  distribution  is  defined  by  Equation  3.15,  the 
extrapolation  procedure  for  soil  is  the  same  as  for  rock  and  concrete 
(Appendix  A).  Thus,  for  2-  or  3-D  motion,  the  function  sin  n  is  re¬ 
placed  by  vn/v  >  and  V  by  v  ,  so  that 


(B.l 6) 


o 


Srp 


♦ 


5-19  v 
s 


ii2_i£  Is. 

»2  v 


As  in  Appendix  A,  it  is  assumed  that  for  a  stress  to  act  on  a  given 
projectile  surface  element  there  must  he  contact  with  the  target  (Sec¬ 
tion  2. It)  and  net  motion  into  the  target  (v  >  0),  otherwise  o  *  0  . 
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Table  B.l.  Typical  S-numbers  for  natural  earth  materials. 


S  _ Materials 

1- 2  Frozen  silt  or  clay,  saturated,  very  hard.  Rock,  veathered, 

low  strength,  fractured.  Sea  or  freshwater  ice  more  than 
10  feet  thick. 

2- 3  Massive  gypsite  deposits.*3  Well-cemented  coarse  sand  and 

gravel.  Caliche,  dry.  Frozen,  moist  silt  or  clay. 

U— 6  Sea  or  freshwater  ice  from  1  to  3  feet  thick.  Medium  dense, 

medium  to  coarse  sand,  no  cementation,  wet  or  dry.  Hard, 
dry,  dense  silt  or  clay.c  Desert  alluvium. 

8-12  Very  loose,  fine  sand,  excluding  topsoil.  Moist,  stiff  clay 
or  silt,  medium  dense,  less  than  about  50  percent  3and. 

10-15  Moist  topsoil,  loose,  with  some  clay  or  silt.  Moist,  medium 
stiff  clay,  medium  dense,  with  some  sand. 

20-30  Loose,  moist  topsoil  with  humus  material,  mostly  sand  and 
silt.  Moist  to  wet  clay,  soft,  low  shear  strength. 

1*0-50  Very  loose,  dry,  sandy  topsoil  (Eglin  AFB).  Saturated,  very 
soft  clay  and  silts,  with  very  low  shear  strengths  and 
high  plasticity. d  Vet  lateritic  clays. 


a 

Taken  from  Reference  3. 

White  Sands  Missile  Range,  Sew  Mexico. 
d  Tonopah  Test  Range,  Nevada,  drj  lake  playas. 

Great  Salt  Lake  Desert  and  bay  mud  at  Skaggs  Island. 
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APPENDIX  C 
NOTATION 


a,  b,  c 
A 
CRH 


Empirical  force  coefficients 

Maximum  cross-sectional  area  of  projectile 

Tangent-ogive  radius  of  curvature,  in  calibers  ("caliber 
radius  head") 


dF 


x’ 


dA 

dF 

z 


g 

I 

k 

L 

M 

n 


Differential  element  of  surface  area  on  projectile 

Differential  force  components  in  the  body-fixed  lateral  and 
axial  (x-  and  z-)  directions,  respectively 

Maximum  projectile  diameter 

Maximum  coarse  aggregate  size 

Function  of  velocity 

Net  force  components  in  the  body-fixed  lateral  and  axial 
(x-  and  z-)  directions,  respectively 

2 

Gravitational  acceleration  (32.2  ft/s  ) 

Transverse  mass  moment  of  inertia  about  CG 
Free-surface  parameter,  r3/rp 
Projectile  length 
Projectile  mass,  W/g 
Exponent 


Nrc  Nose  performance  coefficient  in  rock  and  concrete 

Ng  Nose  performance  coefficient  in  soil 

rc  Cylindrical  radius  of  wake  cavity 

rQ  Cylindrical  radius  of  projectile  at  separation  point 

r  Local  cylindrical  radius  on  projectile  surface  at  a  given 

P  point 

rg  Maximum  perpendicular  distance  from  projectile  axis, 

through  dA  ,  to  target  surface,  for  which  stress  relief 
can  occur  due  to  free-surface  effect 


8l 


Rock  quality  designation 
Soil  penetrability  index 


Absolute  local  velocity  of  any  point  on  projectile  surface 


V  vz 


Outward  normal  component  of  local  velocity 

Local  velocity  components  in  body-fixed  lateral  and  axial 
(x-  and  z-)  directions,  respectively 


Projectile  CG  velocity 

Outward  normal  component  of  CG  velocity 


V  V 
x,  z 


Impact  velocity 

Components  of  CG  velocitv  in  body-fixed  lateral  and  axial 
(x-  and  z-)  directions,  respectively 

Projectile  weight 

Lateral  (body-fixed)  position  in  projectile 
Horizontal  position  in  target 

Components  of  acceleration  in  the  X-  and  2-directions, 
respectively 

Unconfined  compressive  strength  of  intact  rock  or  concrete 
Concrete  strength,  corrected  for  aggregate  size 
Rock  strength,  corrected  for  RQD 

Axial  (body-fixed)  position  in  projectile,  measured  forward 
from  CG 

Vertical  position  (depth)  in  target,  measured  downward 
from  target  surface 

Vertical  position  (depth)  of  CG,  measured  downward  from 
target  surface 

Maximum  (final)  depth  of  penetration  of  nose  tip 


A  dot  over  any  quantity  indicates  differentiation  with  respect  to 
time. 


Qesi 


(Vja\ 


\a0\e 


Cop'/ 
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a  Angle  of  attack 

aQ  Initial  angle  of  attack 

0  Constant,  6.25  *  10^  lb^-s^/ft^ 

C  O 

Y  Constant,  1.5<S  *  10  lb/ft 

6  Lateral  displacement  of  point  on  projectile  axis,  relative 
to  wake  cavity 

Initial  displacement  of  point  on  projectile  axis 

C  Axial  distance  aft  of  the  nose  tip 

CQ  Axial  distance  from  nose  tip  to  separation  point 

CCG  Axial  distance  from  nose  tip  to  CG 

n  Local  angle  between  projectile  axis  and  line  tangent  to 
projectile  surface 

n  Cone  half-angle 

c 

0  Obliquity  angle,  measured  counterclockwise  from  target-fixed 
Z-axis  to  projectile-fixed  z-axis 

0  Angular  velocity 

0  Angular  acceleration 

0q  Initial  obliquity  angle 

u  Constant,  1.96  *  10^  lb/ft 

p  Mas3  density  of  rock  or  concrete 

a  Local  compressive  normal  stress  on  projectile  surface 

a  ,  a  Components  of  local  compressive  normal  stress  in  body-fixed 

lateral  and  axial  (x-  and  z-)  directions,  respectively 

♦  Local  angle  of  approach  between  projectile  surface  element 
and  target,  based  on  CG  velocity  components 

Aftbody  flare  angle 

♦  ,  Angle  of  approach  at  which  separation  occurs  (wake  separa- 

nin  tion  angle) 

Azimuthal  angle  on  projectile 
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